Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



Discovery Potential of the Standard Model Higgs Boson Through 
H WW Decay Mode with the ATLAS Detector at LHC 

H.J. Yang (on behalf of the ATLAS Collaboration) 

Department of Physics, University of l\/licliigan, Ann Arbor, Ml 48109- 1 120, USA 



We report results of a study of the Standard Model Higgs boson discovery potential through the W-pair leptonic 
decay modes with the ATLAS detector at LHC at 14 TeV center-of-mass energy. We used MC samples with full 
detector simulation and reconstruction of the ATLAS experiment to estimate the ATLAS detection sensitivity 
for the reaction of pp — > _ff — > WW — > eufiu with no hard jet or two hard jets in the final states. The prospects 
for the Higgs boson searches at ATLAS are presented, including trigger efficiencies and data-driven methods to 
estimate the backgrounds using control samples in data. With 10 fb~^ of integrated luminosity, one would expect 
to discover a Standard Model Higgs boson with ATLAS detector in Higgs boson mass range 135 < mn < 190 
GeV. If the Higgs boson does not exist, we will extend and confirm the exclusion produced by the Tevatron 
Higgs boson search result, which has ruled out the Higgs boson mass range of 160 < mu < 170 GeV at 95% 
confidence level. If Higgs boson would be discovered, ATLAS could measure its mass with a precision of about 
2 and 7 GeV for Higgs boson mass around 160 GeV and 130 GeV, respectively. 



1. Introduction 

The Standard Model (SM) of particle physics is 
a successful theory confirmed by numerous measure- 
ments since the 1970's. After the direct observations 
of the massive top quark and the Vt at Fermilab in 
1995 and 2000, respectively, only one fundamental 
particle predicted by the SM, the Higgs boson, re- 
mains to be discovered. The role of the Higgs boson is 
critical in the SM because fermions and W/Z bosons, 
via a process of spontaneous symmetry breaking, ac- 
quire their mass through interactions with the Higgs 
boson fields. In the SM the Higgs boson couplings to 
all fermions and electroweak gauge bosons depend on 
their masses. Generically, the SM Higgs boson cou- 
ples most strongly to heavy particles. However, the 
SM doesn't predict the Higgs boson mass, which re- 
sulted in a great experimental challenging to search 
for the Higgs boson in a very wide mass range and 
through many different decay modes. 

The Leading Order (LO) SM Higgs boson produc- 
tion mechanism at the LHC (Large Hadron Collider) 
is shown by the Feynman diagrams in Figure[l] includ- 
ing the processes of the Gluon-Gluon Fusion (GGF), 
the tt fusion, the Higgs-strahlung and the Vector- 
Boson Fusion (VBF). The dominant process is the 
GGF Higgs boson as shown in the next-to-leading- 
order (NLO) cross-section plot in the Figure [2] Higgs 
boson production via VBF is also an important pro- 
duction mechanism for Higgs boson searches at the 
LHC with the specific signature of a Higgs boson with 
two hard jets in the final states. The Higgs boson pro- 
duction cross-section is very small compared to the 
enormous QCD jet productions at the LHC. The abil- 
ity to detect and identify the Higgs boson decay final 
states (see Figure [S]) [Tj is crucial to the success of the 
Higgs boson discovery. 
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Figure 1: SM Higgs boson production LO Feynman dia- 
grams. 



The direct search at LEP sets a lower Higgs boson 
mass limit of 114.4 GeV at 95% confidence level [2]. 
Recent searches at Tevatron exclude the Higgs boson 
mass between 160 GeV and 170 GeV at 95% confi- 
dence level |3J. The electroweak data constraints de- 
duced from consistency conditions of the SM can be 
used to derive an upper limit of SM Higgs boson mass 
of 190 GeV at 95% C.L.[1, 4 . Therefore the most in- 
teresting mass range indicated by direct Higgs boson 
search and indirect constraints is between 114 GeV 
and 190 GeV. 
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Table I Monte Carlo generators and cross sections for 
the Higgs boson signal and background processes. The 
W+jets cross sections listed is the cross section per lepton 
flavor. 



Physics Process 


Generator 


a(pb) 


gg^ WW{Mh = 170 GeV) 


MC@NLO 


19.418 


VBF H WW{Mh = 170 GeV) 


Pythia 


2.853 


qq/qg WW 


MC@NLO 


111.6 


gg WW 


GG2WW 


5.26 


pp — > ti 


MCQNLO 


833.0 


Z ^ TT + jets 


ALPGEN 


2015 


W ^ + jets 


ALPGEN 


20510 



Figure 2: SM Higgs boson production NLO cross-sections 
at the LHC at 14 TeV 
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Figure 3: SM Higgs boson decay branching ratios 



2. Monte Carlo Samples 

As shown in the Figure |2j there are two domi- 
nant Higgs boson production modes: GGF and VBF, 
which are in the kinematic region of interest for the 
H WW decay mode. When both W bosons decay 
to leptons, the Higgs boson signature includes two en- 
ergetic leptons with large missing transverse energy 
in the final state. These two leptons tend to go same 
direction due to spin correlation. For the VBF Higgs 
boson, the final state also includes two hard jets which 
tend to be well-separated in pseudo-rapidity. The 
Higgs boson signal and major background sources (SM 
WW, ti, Z-Hjets and W-|-jets) are hsted in Table |l| 
The MC generators we used to produce MC events 
and corresponding cross sections are also listed in the 
Table H 



3. Leptonic W Pair Production with No 
Hard Jets 

The event selection for H WW ivtv channel 
consists of a set of simple cuts. 

• Require that the selected event has exactly two 
isolated, opposite-charge leptons (electron or 
muon) with Pt > 15 GeV. 

• To suppress single-top background, backgrounds 
from dileptonic decays of bb and cc resonances, 
and lepton pairs from b c cascade decays, re- 
quire that the invariant mass mu of the leptons 
is between 12 GeV and 300 GeV. 

• Require that the event has missing transverse 
energy E^'"' > 30 GeV. 

• To suppress Z tt background, reject events 
which have invariant mass of a hypothetical r 
pair in the range \MrT — Mz\ < 25 GeV. 

• To suppress backgrounds from top quark decays, 
reject events that contain any hard jets with 
PT > 20 GeV and \t]j \ < 4.8. 

• To further suppress the top-related background, 
reject events with any jets with px > 15 GeV 
and a b-tagging weight greater than 4. 

Table [H] shows the cross sections in the e/i channel 
for signal and background after these cuts. 

Events are required to pass at least one of the 
ATLAS single-lepton or double-lepton triggers [1]. 
The following notation is used to label different trig- 
ger items: e (electron), EM (electromagnetic), MU 
or mu (muon), I or i (isolated). As an example 
2EM15I means two isolated electromagnetic objects 
with transverse energy Et > 15 GeV. The Level-1 
(LI) trigger menus used here are 2EM15I, EM25I, 
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Table II Cut fiows(in fb) for Mh = 170 GeV in the H + Oj, Ti" ^ WW ev^jLV channel. The WW background contains 
the two processes qq WW and gg WW . 



Selection 


Selection Cuts 


99 


tt 


WW 


Z ^ TT 


W+jets 




Lepton Selection + 


166.4 


6501 


718.12 


4171 


209.1 


Pre-selection 


EJp"' > 30 GeV 


147.7 


5617 


505.25 


526.3 


181.6 




Z ^ TT Rejection 


145.8 


5215 


485.12 


164.2 


150.4 




Jet Veto 


61.80 


14.84 


238.35 


31.91 


76.12 




b-veto 


61.56 


6.85 


237.87 


30.76 


76.12 


Signal region 


A0 < 1.575, Mt < 600 GeV 


50.6±2.5 


2.3±1.6 


85.4±2.7 


< 1.7 


38±38 


Control region 


A0 > 1.575, Mt < 600 GeV 


10.9±1.1 


4.6±2.3 


151.9±3.6 


30.8±4.2 


38±38 



EM60I, MU20 and MU40. For the Level-2 (L2) trigger 
and Event Filter (EF) , events are required to pass the 
e25i, 2el5i, e60 or mu20i triggers. For the e/i chan- 
nel, the trigger efficiency for LI is 99.0%; for L2 it is 
96.7%, and for EF it is 95.2%. The trigger efficiency 
for signal is quite high and the trigger efficiency does 
not distort the shapes of the kinematic variables of in- 
terest in the signal in a significant way. Table |TT] does 
not include the trigger efficiency. 

After the basic event selection, the signal to back- 
ground ratio is about 1:4. This analysis focuses on 
three variables: 



ATLAS 



Signal 
QCD WW 




• the transverse opening angle A(j)u; a cut on this 
variable exploits differences in the spin corre- 
lations in the WW system in the Higgs boson 
signal and the WW background (Figure |4]) . 

• the transverse momentum of the WW system, 

which tends to be slightly larger for sig- 
nal than for WW background because gluon- 
initiated processes tend to have more initial- 
state radiation than quark-initiated processes 
(Figure |5|. 

• the transverse mass which is defined as M' 



{p^ + pip^''^'^)'^ (Figure 



3.1. Fitting Algorithm 

A 2-dimensional fit of transverse mass (Mt) and 
in two bins of the dilepton opening angle A(l)u 
in the transverse plane has been implemented. After 
the preselection cuts and the additional requirement 
that Mt < 600 GeV is applied, the remaining events 
are separated into two subsamples, one with A0« < 
1.575 (signal region) and the other with A(j>u > 1.575 
(control region), as listed in Table [Tl| 

The top background is estimated with the help of 
b-tagged control samples with the same kinematic 



Figure 4: Transverse opening angle Acpu of two leptons 
after preselection cuts. 
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Figure 5: Transverse momentum p7 of the WW system 
after preselection cuts. 



cuts as the signal-enriched and background-enriched 
regions. The tt cross section is estimated based on 
Monte Carlo samples that use MC@NLO to model 
the top background. Standalone fits are performed 
on the b-tagged control samples before the fit to the 
b- vetoed regions begins, the top background in the b- 
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Figure 6: Transverse mass Mt for events with A(f)u < 
1.575 and > 20 GeV assuming a SM Higgs boson 

with AIh = 170 GeV and 10 fb~^ integrated luminosity. 



vetoed regions is estimated by extrapolating both the 
shape and the normahzation from the b-tagged region 
to the b-vetoed region based on ratios obtained from 
MC@NLO. 

The Z TT background is normahzed and its 
shape is determined by studying a sample of Z ^ fifi 
events taken from real data, where the reconstructed 
muons are replaced by simulated taus. Two muons 
events with a dimuon invariant mass between 82 and 
98 GeV are selected, and the same jet veto is ap- 
plied to selected events. After event selection cuts are 
applied, the effective cross section of Z ^ tt back- 
ground is about 250 fb. A standalone fit to these 
"data-Monte-Carlo" events is performed to determine 
the shape and normalization of the Z —^ tt back- 
ground; in the final fit, the shape parameters obtained 
from the fit to control sample are fixed and the nor- 
malization is rescaled by a factor that is assumed in 
this study to be well-predicted. 

W-fjets is one of the main sources of fake back- 
grounds for the dilepton channels, it is crucial to 
achieve a good rejection against this background. The 
average jet fake as electron rate is (1.7 ± 0.2) x 10"'' 
before isolation and (6.7 ± 1.5) x lO"'^ after electron 
isolation. The corresponding jet fake as muon rate 
is (1.7 ± 0.6) X 10~^ after isolation cuts are applied. 
For W-|-jets background estimation, due the limited 
size of the available W-l-jets Monte Carlo sample, the 
shape of the transverse mass and pjf^ distributions 
are taken from a set of events with loosen isolation 
and shower shape cuts. Further detailed studies are 
needed to estimate W-|-jets background. 

Once the fits to the control samples are completed, 
a simultaneous fit to the two A0££ bins in the b-vetoed 
region is performed. A few of the parameters that de- 
scribed the shape of the transverse mass and pj^^ 



distributions of the WW background are allowed to 
float in the fit. The normalizations of the WW back- 
ground are free to float independently. However, we 
add a penalty term of the form (Rfu — Rtrue) / '^w 
where Rfu is the ratio of the best-fit number of WW 
background events in the small-Ac^^^ region over the 
number in the large large-A(/)^£ region, Rtme is the 
Monte Carlo prediction of the ratio taken from the 
central- value calculation, and an is the uncertainty in 
the prediction of Rtme, taken to be 10%. 
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Figure 7: The log Likelihood Ratio distributions for 
background-only Monte Carfo in H -I- Oj, H ^ WW 
eu^v assuming 10 fb~^. 
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Figure 8: The pull distributions for Mh = 170 GeV for 
Higgs boson signal and background Monte Carlo in H -|- 
Oj, H — > WW ev^Jiv assuming 10 fb~^. 



In order to demonstrate the robustness of the fit 
against systematic uncertainties, top Monte Carlo has 
been used to compute the sampling distributions for 
the Likelihood Ratio in several scenarios where the 
"true" probability distribution has been distorted to 
model various sources of systematic error. Seven 
distorted scenarios are considered: four altered 
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scale choices (factorization and renornialization scales 
raised and lowered by factors of 8), two alternative 
top background models (based on leading-order pp —^ 
WWbb and leading-order pp ^ tt ^ WWbb), and 
one alternative model of all irreducible backgrounds 
where the x and y components of E'™*^* have been in- 
dependently smeared by 5 GeV each. Figure [7] shows 
the Likelihood Ratio distributions for background- 
only outcomes and Figure [8] represents the distribu- 
tions of pulls of the fitted Higgs boson mass for signal 
plus background with a true Higgs boson mass of 170 
GeV. 

The linearity of the mass determination as a func- 
tion of the true Higgs boson mass is shown in Figure [9j 
The line shows the mean of a Gaussian fit to the re- 
gion around the peak of the distribution of best-fit 
Higgs boson masses in the Monte Carlo sample for 
the case of nominal detector performance. The green 
band shows the width of the Gaussian fit and is a di- 
rect measure of the variability of the mass estimate on 
repetition of the experiment, the error bars show the 
median fit error. The typical variability of the mass 
determinations with 10 fb~^ integrated luminosity are 
5.2 GeV, 1.6 GeV and 4.2 GeV at Mh = 130, 160 and 
190 GeV, respectively. 
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Figure 9: The linearity of the Higgs boson mass determi- 
nation for H -I- Oj, -ff ^ WW eufiu at 10 fb"^ inte- 
grated luminosity. 



The H WW eiyfiv is the most promising chan- 
nel for Higgs boson masses near the WW threshold. 
The expected SM Higgs boson detection significance 
assuming 10 fb"-'^ integrated luminosity is shown in 
Figure [TO] 
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Figure 10: The expected Higgs boson detection signifi- 
cance for the H -I- Oj, ii" ^ WW eufif at 10 fb"^ 
integrated luminosity. 



4. Leptonic W Pair Production with Two 
Hard Jets 

For SM Higgs boson produced through VBF pro- 
duction mechanism, the final state has two energetic 
leptons from both W bosons leptonic decays associ- 
ated with two hard jets. The distinctive characteris- 
tics of VBF Higgs boson signal include: 

• The two jets arising from struck quarks tend to 
be the highest-pT jets in the events, and they 
tend to be well-separated in pseudo-rapidity; 

• they tend to have a large invariant mass; 

• there is very little jet activity in the region be- 
tween the two tagging jets. 

Figure [TT] shows the distribution of the pseudo- 
rapidity gap between two tagging jets (top left), in- 
variant mass of two tagging jets (top right), azimuthal 
angle between two tagging jets (bottom left) and 
transverse energy of the third jet (bottom right). 
Where black solid histograms represent VBF Higgs 
boson signal with Mh = 170 GeV, green dash-dotted 
histograms mean WW background, red dotted his- 
tograms are tt and blue dashed histograms are Z-|-jets 
background. 

The following cuts have been applied to selected H 
+ 2j, H ^ WW Mv events. 

• Two isolated leptons with px > 15 GeV 

• Missing transverse energy i?™*'"' > 30 GeV 

• At least two jets withpT- > 20 GeV and \ri\ < 4.8 

• The two jets with highest transverse momentum 
are required to be in opposite hemispheres with 
A77(jeil,jei2) >3 
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Figure 11: Pseudo-rapidity gap between two tagging jets(top left), invariant mass distribution of two tagging jets(top 
right), azimuthal angle gap between two tagging jets(bottom left) and Et of the third jet(bottom right). 



Table III Cut flows(in fb) for Mh = 170 GeV in the H + 2j, -> WW ^ eivfiu channel. 



Selection Cuts 


VBF H (170 GeV) 


tt 


WW-l-jets 


Z ^ TT 


W-l-jets 


Lepton Selection 


30.20 


8317 


838.96 


2096 


1323 


Forward Jet Tagging 


17.27 


946.6 


32.77 


79.3 


31.83 


Leptons Between Jets 


16.47 


617.8 


22.92 


55.13 


27.91 


Z ^ TT Rejection 


15.68 


561.8 


21.20 


39.03 


27.91 




12.78 


425.9 


15.28 





13.96 


b-veto 


12.67 


206.72 








signal box, b-jet veto 


9.28±0.27 


28.5±5.7 


4.75±0.30 




4.3±4.3 


signal box, no b-jet veto 


9.65 


114.2 


4.99 




6.07 


Control box, b-jet veto 


3.02±0.15 


89±10 


9.78±0.43 




7.9±5.0 


Control box, no b-jet veto 


3.13 


311.7 


10.28 




7.89 



• Require that both leptons are between the two 
leading jets in pseudo-rapidity 

• Reject event which has \MrT — Mz\ < 25 GeV 

• Require transverse mass 50 < Mt < 600 GeV 

• Require transverse mass m^" > 30 GeV, m^^ = 

\/2Pt(«)^'t • (1 - cosA(j)) where is the 
angle between the di-lepton vector and the 
pmzss yQQf^Qj- the transverse plane. 



• To suppress b-related background, b-veto cuts 
are applied on both tagging jets 



Table [ml shows the cut flows (in fb) for Mh = 170 
GeV in the H + 2j, ^ WW ev^iv channel. If 
the event has A(/)« < 1.5 and At;^ < 1.4, it lies in 
the signal box; otherwise, it lies in the control region. 
The trigger efliciency of 99.0% after Level 1, 96.8% 
after Level 2, and 94.5% afte the Event Filter. 

After the preselection, a four- variable Neural Net- 
work is used to further enhance the separation be- 
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Figure 12: The Neural Network output distribution in the 
signal box, for the 170 GeV H + 2j, // ^ WW ^ eufiu 
at 10 fb~^ integrated luminosity. 
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Figure 14: The expected Higgs boson detection signifi- 
cance for the H + 2j, H ^ WW eufif at 10 fb"^ 
integrated luminosity. 
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Figure 13: The transverse mass distribution for events in 
the signal box with Neural Network output larger than 0.8 
for the 170 GeV R + 2j, H ^ WW -> eiyfiu at 10 fb^^ 
integrated luminosity. 



curve represents 2-lepton background and green curve 
is faked background from W+jets. 

The fit is a two-dimensional fit to the Neural Net- 
work output and transverse mass distributions. Both 
the signal model and the background model are uncor- 
related product probability density functions (PDFs) . 
The Neural Network output distribution for the back- 
ground in the signal box is taken to be the same as the 
distribution in the control region, but it is multiplied 
by a linear extrapolation factor. Apart from the slope 
of this extrapolation factor, all parameters governing 
the shape of the Neural Network output distribution 
in the two regions are required to be the same. The 
expected significance is shown as a function of the true 



Higgs boson mass in Figure 14 



tween the signal and the background. The inputs to 
the Neural Network are: 

• Ar](jetl, jet2), the pseudo-rapidity gap between 
two tagging jets 

• Mjj , the invariant mass of two tagging jets 

• the transverse momentum of the leading non- 
tagging jet in the region I77I < 3.2 

m rj* = 773— (77i-|-?72)/2, the pseudo-rapidity gap be- 
tween two tagging jets and the third non-tagging 
jet. 

The Neural Network output distribution in the sig- 



nal box is shown in Figure 12 The transverse mass 



distribution for events in the signal box with Neu- 
ral Network output larger than 0.8 is shown in Fig- 
ure |13[ where black dots with error bar show Pseudo- 
data, blue curve shows best-fit Higgs boson signal, red 



5. Further Improvements 



We expect further improvements for Higgs boson 
detection sensitivity. The analysis presented above 
only includes the ev^LV channel with jet-veto or two 
forward jet tagging in events. Including the evev and 
jivjiv channels and adding dilepton events with 1 jet 
in final states will certainly increase the detection sen- 
sitivity in ATLAS. Based on LEP and Tevatron Higgs 
boson search experiences, the analysis can be carried 
out by using mativariate techniques, such as Artificial 
Neural Networks (ANN) and Boosted Decision Trees 
(BDT) which would further improve the detection 
sensitivities. ATLAS physics group has extensively 
explored these advanced analysis methods [1] and fully 
developed the tools that will be used in Higgs boson 
search program with LHC collision data. 
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Figure 15: The expected Higgs boson detection signifi- 
cance for tfie H + Oj and H + 2j, _ff — > WW — > ei/nv at 
10 fb~^ integrated luminosity. 



grated luminosity, we expect to reach 5cr discovery if 
SM Higgs boson mass between 150 GeV and 180 GeV. 
The combined mass determination for the combined 
fit of H + 0/2j is shown in Figure 16 The correspond- 



ing combined significance as a function of Higgs boson 



mass is shown in Figure 15 
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Figure 16: The linearity of the mass determination for the 
combined fit of H + 0/2j, H WW eufiu at 10 fb"^ 
integrated luminosity. 



6. Conclusions 

The prospects for SM Higgs boson searches in the 
WW decay mode have been studied using a real- 
istic model of the ATLAS detector. The H -I- Oj, 
H — > WW — > evfiv channel is very promising for 
Higgs boson masses in the region around the WW 
threshold. With 10 fb~^ of integrated luminosity, one 
would expect to be able to achieve a 5a discovery with 
the H WW — *■ ei'fii' channel alone if there a SM 
Higgs boson with a mass between ^140 GeV and ^ 
185 GeV. A measurement of the mass of the Higgs 
boson at 10 fb~^ integrated luminosity would have 
a precision of less than 2 GeV for Higgs mass bo- 
son of 160 GeV, or a precision of less than 4 GeV 
for Higgs boson with mass of 140 GeV. The H -|- 
2j, if — > WW — > eviiv channel has a smaller event 
rate than the H -I- Oj channel. With 10 fb^^ of inte- 
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